0 5 


偏振 


偏振 
检 侦 


基于 连续 旋转 波 片 调制 的 传 里 时 解 调 方法 研究 


WAF’, RAE 
(1. 中 国 科学 院 云 南天 文 台 ， 云南 昆明 650216; 2 . 中 国 科学 院 大 学 ,北京 100049) 


摘要 : 本 文 给 出 了 在 连续 式 旋转 波 片 调制 模式 下 的 傅 里 叶 分 析 的 推导 过 程 和 解 调 公 
Ky, 并 利用 理论 偏振 轮廓 对 公式 的 准确 性 进行 了 验证 。 分 别 采 用 了 傅 里 叶 分 析 和 解 调 矩 
阵 两 种 解 调 方法 定量 模拟 了 在 连续 式 调制 模式 下 波 片 的 初始 方位 角 误差 ,旋转 定位 误差 
以 及 探测 器 曝光 过 程 中 时 间 差 占 比 三 个 因素 对 测量 精度 的 影响 。 主 要 结论 是 : (1) 简单 
的 倩 里 叶 分 析 不 再 适用 于 连续 式 调 制 的 结果 解 调 .文中 给 出 的 新 倩 里 叶 分 析 方 法 与 解 调 
珑 阵 方法 在 进行 误差 分 析 时 ， 可 以 得 到 基本 一 致 的 结果 。 (2) 考虑 初始 方位 角 误差 引起 
的 影响 ,我 们 发 现 对 于 线 偏 振 信 号 来 说 ， 采 用 傅 里 叶 分 析 和 解 调和 矩阵 可 得 到 非常 一 致 的 
结果 ， 而 对 于 圆 偏振 来 说 ， 基 于 解 调和 矩阵 的 方法 得 到 的 结果 误差 相对 较 小 。 但 是 两 种 解 
调 方法 都 反映 出 初始 方位 角 误 差 对 圆 偏振 和 线 偏振 产生 的 影响 是 一 致 , 且 相 对 误差 的 大 
小 与 偏振 信号 本 身 的 强 弱 有 关 。 初始 方位 角 的 误差 要 求 在 十 几 角 秒 量 级 才能 满足 结果 的 
相对 误差 在 10 SBF. (3) 考虑 旋转 定位 误差 引起 的 影响 时 , MAMA +A 
析 和 解 调 矩 阵 得 到 的 结果 非常 接近 。 两 者 同时 都 显示 出 旋转 定位 误差 对 线 偏振 信号 造成 
的 影响 更 为 明显 ， 且 相对 误差 的 大 小 也 与 信号 的 强 弱 有 关 。 当 波 片 的 重复 定位 精度 在 
10" 左 右 时 ,对 10-? 量 级 的 弱 偏 振 信 号 来 说 ,测量 误差 也 能 在 10-3 的 量 级 左右 。 从 这 一 点 
上 ， 我 们 发 现 连续 式 调制 模式 对 波 片 旋转 的 定位 精度 要 求 要 明显 高 于 步 进 式 调制 模式 。 
(4) 两 种 解 调 方法 的 结果 均 显示 ， 在 连续 式 调制 模式 下 ， 波 片 的 调制 周期 与 探测 器 曝光 
时 长 的 时 间 差 (A) 会 造成 线 偏振 信号 之 间 的 串扰 ， 即 时 间 差 占 比 (A/T) 对 线 偏振 信号 造 
成 的 影响 比 对 圆 偏振 信号 的 影响 更 为 明显 。 当 A/T 小 于 1% 时 ， 线 偏振 信号 相对 误差 达 
到 10-3 量 级 。 
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偏振 辐射 测量 是 获得 天 体 物理 状态 的 一 个 重要 手段 , 光源 或 传播 路 径 的 各 向 异性 是 导 
产生 的 主要 机 制 , 因此 偏振 辐射 可 以 提供 有 关 诸 如 天 体 人 磁场 或 辐射 场 等 矢量 场 的 重要 信 


Bo 在 太阳 实测 物理 中 , 利用 某 些 原子 谱 线 进行 高 精度 偏振 辐射 测量 是 目前 探测 太阳 矢量 磁 
场 的 一 种 主要 手段 。 偏 振 测量 精度 通常 要 求 在 10 或 更 高 的 量 级 。 


我 们 通常 采用 斯 托 克 斯 偏振 矢量 ， 即 Sm = U, Q, U,V)" (TREE) 来 描述 一 束 光 的 偏振 
。 因 为 探测 器 只 能 测量 强度 信号 ， 因 此 我 们 需要 将 4 个 偏振 矢量 进行 线性 组 合并 进行 多 


测量 (或 称 之 为 偏振 调制 )， 最 后 才能 从 测量 的 强度 信号 中 反 推 出 完整 的 偏振 信息 (或 称 为 


解 调 ) 。 实 现 偏振 矢量 线性 组 合 的 光学 元 件 组 称 之 为 偏振 分 析 器 ， 通 常 由 偏振 调制 器 和 
器 两 部 分 组 成 。 偏 振 调制 过 程 的 数学 表达 可 表示 为 as = 0 Sin。 其 中 0 是 一 个 n x 4 的 
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和 矩阵 (n 为 测量 次 数 ，n > 4), 被 称 为 调制 矩阵 。 解 调 过 程 则 是 调制 过 程 的 逆 过 程 ， 可 表达 为 
D> lops = Sam。 这 里 D 称 为 解 调 矩阵 ， 是 调制 矩阵 的 逆 ”。 利 用 正确 的 调制 及 解 调和 矩阵 ， 我 
们 就 可 以 根据 观测 计算 出 输入 的 偏振 信和 号 


波 片 或 波 片 组 合 是 一 种 常见 的 偏振 调制 器 元 件 , 利 用 波 片 来 实现 偏振 调制 也 有 多 种 形式 ， 
例如 , 在 电机 带动 下 旋转 固定 延迟 量 的 光学 波 片 的 快 轴 ， 或 者 通过 改变 电压 来 改变 液晶 波 片 
的 延迟 量 或 光 轴 方向 等 等 。 相 对 使 用 液晶 波 片 来 说 ,使 用 光学 波 片 的 优点 主要 体现 在 以 下 几 
A: (1) 性 能 比较 稳定 ， 在 有 效 通 光 孔径 范围 内 ， 波 片 的 面 形 精度 以 及 延迟 量 分 布 均匀 性 
相对 较 好 。(2) 波 片 的 延迟 量 可 以 根据 科学 需求 进行 优化 。 优化 原则 通常 是 考虑 到 各 种 科学 
目标 对 不 同 偏振 矢量 的 测量 效率 要 求 的 不 同 。(3) 适 用 于 较 宽 的 波长 范围 ， 消 色差 波 片 可 以 
实现 多 波段 的 同时 观测 ,因此 多 个 太阳 望远镜 采用 了 可 旋转 光学 波 片 作为 偏振 调制 器 ,例如 ， 
地 面 设备 有 美国 Vacuum Tower Telescope 的 Advanced Stokes Polarimeter (ASP) All 1.6 
米 口 径 的 Goode Solar Telescope (GST) “， 中 国 的 一 米 新 真空 太阳 望远镜 (NVST)"” 和 光 
纤 阵 列 太阳 望远镜 (FAS0T) “等 等 。 空间 设备 有 日 本 的 Hinode 卫星 的 SOT 仪器 搭载 的 SP 载 
荷 和 探 室 火箭 CLASP2” 等 。 用 于 日 冕 磁场 测量 的 印度 Visible Emission Line 
Coronagraph (VELC) “日 园 仪 也 是 采用 旋转 波 片 作为 偏振 调制 器 。 从 上 述 仪器 的 偏振 调制 方 
式 中 我 们 发 现 , 波 片 旋转 可 采用 两 种 模式 ,一 种 是 波 片 步 进 旋转 ,探测 器 在 波 片 每 步 旋转 到 
位 并 停止 时 进行 数据 采集 (例如 NVST 和 FAS0T)。 另 一 种 是 波 片 连续 旋转 ,探测 器 同步 采集 
波 片 旋转 过 程 中 的 强度 变化 (例如 GST, Hinode/SP, VELC 等 ) 。 

在 两 种 模式 中 ， 等 间隔 采集 是 一 种 常见 方式 ， 即 波 片 的 旋转 角度 为 86 = 27/1, n 是 在 波 


片 旋转 一 周 所 需要 的 步 数 或 探测 器 采集 的 次 数 。 以 步 进 式 情况 为 例 , 探测 器 采集 到 的 强度 信 
号 Joos 与 波 片 旋转 角 9、 波 片 延迟 量 6 以 及 入 射 偏振 信号 Sm 之 间 的 关系 为 ”: 


1 
1ops(9) = 5 [( + aQ) + bQcos46 + bUsin46 + cVsin26]#(1) 


1 + cosé 1 — cos 
a= 5 ,b= 5 > c=-—sind; #(2) 
我 们 发 现 ，Q@、V 以 相同 的 频率 (4 倍 波 片 旋转 的 频率 ) 被 调制 进入 fons ， 但 两 者 相位 相差 
= TT/2。V 则 以 2 倍 波 片 旋转 的 频率 调制 进入 10ps。Q、U 和 V 的 调制 振幅 ( 即 a、b 和 c) 也 仅 与 波 片 
延迟 量 65 有 关 。 因 此 ,除了 使 用 解 调 矩阵 ( 即 利用 前 述 公 式 Sin = D1ops) 计 算出 入 射 的 Si 外 ， 
我 们 还 可 以 根据 傅 里 叶 分 析 解 调 出 $i,。 即 ， 


2 
Oa 2d SORE {FFT [tovs (f = Sp /b 
2 
U = ims jh = 2m {rrr | (f = Si /b #(3) 


V = 2Im2/c = 2Im{FFT[Ion,f = pure 


这 里 符号 Re 和 1m 分 别 代 表 傅 里 叶 分 析 频 率 的 实 部 和 虚 部 ， 而 数字 ‘4? 或 者 ‘2* 则 代表 对 应 
的 傅 里 叶 分 析 频 率 ， 并 与 f 的 值 相对 应 。 


但 是 在 连续 式 旋转 模式 下 , 公式 (1) 需要 进行 调整 , 即 进行 角度 积分 ( 详 见 下 文 的 公式 推 
导 ) ， 此 时 上 文 提 到 的 简单 的 傅 里 时 解 调 公式 就 不 再 适用 了 。 因 此 针对 连续 式 调 制 ， 人 们 还 
是 普遍 采用 基于 解 调 矩 阵 D 的 解 调 方法 “””。 


相对 于 步 进 式 来 说 , 连续 式 波 片 旋转 在 调制 频率 上 具有 明显 的 优势 , 在 国际 上 得 到 普遍 
使 用 。 因 此 ,本 文 的 研究 工作 主要 是 针对 连续 式 波 片 旋转 模式 展开 ， 首 先 推导 出 了 该 调制 模 


E 


= 


式 下 精确 的 傅 里 叶 分 析 的 解 调 公式 , 并 利用 理论 模拟 验证 了 公式 的 正确 性 。 随后 利用 两 种 不 
同 的 解 调 方法 ( 解 调和 矩阵 以 及 传 里 叶 分 析 ) ， pees 波 片 的 初始 方位 角 误差 、 ede abana 
差 以 及 探测 器 采集 过 程 中 的 时 间 差 占 比 等 因素 对 解 调 结果 的 影响 。 最 后 给 出 讨论 分 析 和 结论 


1 公式 推导 与 验证 


在 本 节 中 , 我 们 首先 给 出 连续 式 旋转 波 片 的 傅 里 叶 解 调 公 式 的 推导 过 程 , 然后 理论 模拟 
已 知 偏振 矢量 的 连续 式 调制 观测 , 之 后 分 别 利 用 新 推导 出 的 公式 以 及 步 进 式 情况 下 的 公式 对 
模拟 观测 结果 进行 傅 里 叶 分 析 , 最 后 将 解 调 结果 与 已 知 的 输入 信号 进行 比较 , 验证 公式 的 正 
确 性 。 
1.1 公式 推导 


在 连续 式 模式 下 ， 探 测 器 单 次 曝光 采集 得 到 的 强度 为 : 


Il, = OF 1+ OQ + OLU+ OV #(4) 
HF, OWENE, HEREN: 


;_ 1% 
Oj "arh =! 


OS al, (a + b cos40)d0 #(5) 
了 三 
07 = To al, b sin46d6é 


天 二 1 | 
Or = T c sin26 dé 


其 中 aq、b、c 的 定义 与 公式 (2) 相同 。j = 1，2...n (n 代表 在 波 片 旋转 一 周 内 探测 器 的 
采集 总 次 数 ) 。 测 器 第 /次 曝光 过 程 中 ， 波 片 快 轴 从 67 转 到 9 ， 期 间 波 片 快 轴 扫 过 的 角度 
为 Ag7 = 6; 一 91 。 当 探测 器 曝光 采集 时 间 相 等 ， 即 均匀 采集 时 ，467 = A9 = 2x/n。 如 果 
假设 波 片 JALE 与 检 偏 器 的 透 过 轴 方 向 存在 一 定 的 夹 角 & ,那么 了 = a + G- 1)40. 
将 积分 公式 计算 后 带 入 公式 (4) ， 可 以 得 到 观测 到 的 强度 轮廓 为 以 下 表达 式 : 


I, =1+ aQ 


1 P Ge + w „ bsin(2A6/) sin(240/ + 4a) 


+ cos(461) <5 5 5 
了 J š ive) F 
+sin(46i) 5 1 pcan na +), ee 
— sin(20/ 1) = sin(ao/) cos(A@/ + 2a)V 一 cos(26)) — sin(A0/) sin(Ad/ + 2a)V#(6) 


Ags 

分 析 上 述 公式 我 们 发 现 ，@，U,，Y 信 号 不 再 如 前 文 的 公式 0) - (3) 所 示 的 那样 与 菜 种 调 
制 频率 存在 一 个 简单 的 对 应 关系 , 例如 cos46; 的 系数 中 既 包 含 了 @ 信 号 也 包含 了 U 信 号 。『 信 
号 仍然 只 与 cos291 的 系数 有 关 。 因 此 ， 经 过 一 系列 的 代数 运算 ， 我 们 推导 出 了 新 的 传 里 叶 
分 析 的 解 调 公式 如 下 : 


„ Cos(2A07) Re4 一 sin(2A6/ + 4a) Im4 
~ bsin(2A0/) 
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,Sin(2A07 + 4a)Re4 + (2A67 + 4a)Im4 


b sin(2A07) 
Im2 


其 中 ，Re4、Im4 和 Im2 的 含义 与 公式 (3) 中 的 一 样 。 
1. 2 公式 准确 性 的 仿真 验证 
接 下 来 ， 我 们 利用 理论 模拟 信号 对 公式 (7) 的 解 调 方法 的 准确 性 进行 验证 ， 思 路 如 下 : 


(1) 利用 Rybicki-Hummer (RH) 偏振 辐射 转移 程序 ”模拟 生成 理论 Fel 原 子 谱 线 的 
偏振 轮廓 。 计 算 时 采用 一 维 平面 平行 大 气 假设 ， 宁 静 区 FALC 大 气 模型 以 及 典型 磁场 位 形 ( 磁 
场 强 度 1000 高 斯 ， 倾 斜 角 45°*， 方 位 角 45°) 。 工 作 谱 线 FeI 中 心 波长 630. 15 nm。 理 论 计算 偏 
振 轮 廓 如 图 1 中 实 线 所 示 。 这 里 需要 指出 的 是 ， 我 们 使 用 偏振 轮廓 Sin (A) 而 不 使 用 一 组 特定 
的 偏振 矢量 5;,， 因 为 不 同 的 波长 点 具有 强 、 弱 不 同 的 偏振 信和 号， 可 以 更 好 的 验证 解 调 公式 
的 普 适 性 。 

(2) 模拟 理想 的 8 步 连续 式 波 片 调制 模式 。 即 在 波 片 旋转 半 圈 ( 即 一 个 调制 周期 ) 的 过 程 
中 探测 器 均匀 采集 8 次 。 此 时 A6 = 22.5 ， 波 片 延迟 量 选择 6 = 127 。 根 据 公式 (6) 计算 所 得 
到 波 片 旋转 1 圈 过 程 中 的 强度 轮廓 is 变化。 分别 利 用 公式 (3) 和 (7) 对 理论 强度 轮廓 jms 进行 
傅 里 叶 分 析 方 法 来 解 调 。 


(3) 对 比 已 知 偏振 轮 廊 与 两 种 解 调 结果 。 


= A 
wae: 


aa 
wavelenginA wating A 


图 1 ie FE ES f ATE TE MESE, MEST, opoillie IIe Q/1. U/I V/I 轮 
Ro SAG CE ite FE EB, AAAH ZS SK (3) 中 简单 傅 里 叶 解 调 得 到 的 结果 。 左 上 子 图 中 ， 虚 线 代 表 傅 里 
叶 解 调 的 强度 轮廓 与 理论 强度 轮 廊 的 差 值 。 


Fig. 1 Comparison between the synthetic Stokes profiles (solid line) and the Fourier demodulation 


results (dots). Here we demodulate Stokes profiles by using Eq. (3). From left to right and top 
to bottom, we display the profiles of I/Ic, Q/I, U/I, V/I, respectively. In the left-top panel, 
the dashed line (red) illustrates the residual between the synthetic Stokes I/Ic profile and the 


demodulated one. 


首先 理论 轮廓 显示 出 , 在 我 们 特殊 的 磁场 位 形 的 假设 下 , CQ/ 信和 号 最 弱 ( 最 大 值 约 0.02) , 


J \ 
aa 


>- #(7) 
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UYI 稍 强 ( 大 约 是 QI 的 3 倍 )，V/YI 最 强 ( 最 大 值 约 0.15)〉。 


其 次 ,我们 在 图 1 中 先 展示 了 采用 简单 傅 里 叶 分 析 ( 即 公式 (3) ) 得 到 的 结果 与 已 知 偏振 轮 
廓 的 对 比 。 各 子 图 分 别 展 示 了 1/Ic、Q/I、U/YI 和 VYI 的 情况 。 实 线 为 输入 的 理论 偏振 轮 廊 ， 
点 为 傅 里 叶 解 调 结果 。 对 比 输入 值 和 结果 我 们 明显 发 现 : 在 理想 的 情况 下 ,简单 的 傅 里 叶 分 
析 的 解 调 结果 与 输入 的 理论 轮廓 明显 不 同 ， 特 别 是 在 线 蕊 以 及 近 线 器 (AAx0.05 nm) 的 强 偏 
振 信 号 区 域 。 差 异 最 明显 的 情况 发 生 在 最 弱 的 @A 轮廓 中 ， 甚 至 发 生 了 符号 反 转 。 而 LT、 
V/I 轮廓 的 最 大 差异 的 相对 值 也 到 25%。7A 轮廓 看 似 相同 ， 但 从 两 者 的 残 差 上 看 出 他 们 最 
大 相对 差异 也 有 接近 2% 。 
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图 2 在 不 同 的 采集 帧 数 情 况 下 ， 已 知 偏振 轮廓 与 简单 傅 里 叶 解 调 轮廓 的 对 比 。 左 CAT) E: 0/1 (V/1) 轮 
廓 对 比 ， 采 集 帧 数 分 别 为 8、32 和 128。 下 : n=128 时 ， 已 知 轮廓 与 解 调 轮廓 的 残 差 。 


Fig. 2 Comparison between the synthetic Stokes profiles (solid line) and the Fourier demodulation 


results when the number of samples per half-rotation of the waveplate increases. Left (right) -upper 
panel: The Q/I (V/I) profile. Bottom: The residual of the synthetic profile and the demodulated 


results with the number of samples equal to 128. 


为 了 分 析 上 述 误差 与 一 个 周期 内 的 采集 次 数 的 关系 ， 我 们 将 采集 次 数 从 8 次 增加 到 32 和 
128。 结 果 如 图 2 所 示 。 图 2 左 、 右 子 图 分 别 展示 了 QV/I 信 号 和 VY/I 信 和 号。 我 们 发 现 随 着 采集 帧 
数 的 逐渐 增加 ， 输 入 和 输出 信号 的 残 差 有 所 减 小 。 当 N=128 时 ，Q/I 信 号 的 残 差 最 大 依然 能 
到 3X10 (是 偏振 信号 本 身 的 10%)，V/I 信 号 的 残 差 非常 小 ， 最 大 约 为 7. 5X10 ，( 是 偏振 信 
号 本 身 的 0.04%) ,但 非常 明显 的 是 , 两 者 残 差 的 轮廓 都 与 偏振 信号 本 身 的 轮廓 非常 相似 ,并 
不 是 随机 误差 。 

接 下 来 ,我 们 采用 新 的 傅 里 叶 分 析 方 法 (公式 7) ) 进行 解 调 ， 结 果 展 示 在 图 3 中 。 解 
调 结 果 与 理论 输入 值 的 差异 与 图 1 构成 了 鲜明 的 对 比 。 为 了 突出 差异 ， 我 们 在 每 张 子 图 中 给 
出 了 两 者 的 残 差 。 从 残 差 来 看 (虚线 )， 解 调 后 的 轮廓 与 输入 轮廓 几乎 没有 明显 差异 ,在 整个 
轮廓 范围 内 残 差 都 维持 在 一 个 极 低 的 误差 范围 (10-16~10-17) 。 
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图 3 己 知 偏振 轮廓 与 新 傅 里 叶 解 调 轮廓 的 对 比 。 实 线 为 已 知 输 入 偏振 轮廓 ， 点 为 利用 公式 (7) 得 到 的 新 傅 
里 叶 解 调 结果 ， 点 - 短 实 线 (红色 ) 为 两 者 的 残 差 。 
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Fig. 3 Comparison between the synthetic Stokes profiles (so 
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results(dots). Here the demodulated Stokes profiles are retrieved using the formulas 


The dash-dotted li 


ne shows the residual between these two profiles. 


RA 


-> 


in Eq. (7). 


图 1、2 的 结果 与 图 3 中 的 结果 对 比 充分 展示 了 简单 的 步 进 式 傅 里 叶 分 析 公 式 已 经 不 再 适 


用 于 连续 式 波 片 旋转 调制 的 解 调 , 而 本 文 


调 结果 。 因 此 ， 下 文中 所 提 及 


HA 


E 导 出 的 新 傅 里 叶 解 调 公式 则 能 非常 准 而 


的 傅 里 叶 分 析 或 个 


EE 叶 解 调 均 是 指 的 公式 (7) 中 给 


2 不 同 解 调 方法 与 连续 式 调制 过 程 中 的 误差 分 析 


从 公式 (6) 可 以 看 出 ， 有 几 个 变量 的 误差 会 影响 连续 式 旋 转 波 片 的 调制 结果 ， 


的 初始 方位 q 角 ( 波 片 的 快 轴 方 向 与 检 偏 器 的 透 过 轴 方 向 的 夹 角 )、 旋 转 波 片 的 定 


得 给 出 解 


的 方法 。 


包括 波 片 


MRE. BK 


片 旋转 速度 与 探测 器 曝光 时 间 不 匹配 等 等 。 本 节 中 , RATASE PE id EE EL op 
种 不 同 的 解 调 方法 ， 定 量 分 析 上 述 变量 的 误差 对 解 调 结果 的 影响 。 


2.1 初始 方位 角 误 差 


光路 中 ， 波 片 的 快 轴 与 检 人 4 


位 角 误 差 。 


误差 的 影响 时 ， 我 们 的 分 析 思 


路 如 下 : 


局 器 的 透 过 轴 方 向 会 存在 一 定 的 夹 角 c， 称 之 为 波 片 的 初始 方 
在 实际 工作 中 ， 我 们 无 法 知道 a 的 精确 值 ， 因 此 产生 一 定 的 计算 误差 。 在 模拟 该 


C) 输入 已 知 偏振 信号 ，5j = (1,Q/1,U/LV/1) = (1,0.02,0.06,0.2)", 分 别 代表 弱 、 
中 、 强 偏振 信号 (此 后 QJ/1、U/1、V/I 信 号 统一 用 Q、U、V 信 号 指 代 ) 。 模 拟 8 步 连 续 式 调制 
过 程 , 设 波 片 初始 方位 角 为 a, 在 探测 器 单 次 曝光 过 程 中 波 片 转动 角度 r/8 , 波 片 延迟 量 127°。 
根据 公式 (6) 生 成 所 观测 到 的 强度 轮廓 10ps。 


(2) 采用 解 调 矩阵 和 侍 里 叶 分 析 方 法 分 别 对 观测 到 的 强度 轮廓 进行 解 调 ， 对 比 解 调 结 


果 与 输入 结果 的 相对 误差 ， 定 义 为 |5iw 一 Sout |/Sine 


内 为 a 值 大 小 未 知 ， 所 以 计算 解 调 矩阵 


或 傅 里 叶 分 析 时 都 设 为 a = 0? 的 理想 情况 。8 步 连续 式 调制 的 理想 调制 矩阵 0 如 下 ,其 解 调 矩 


公式 (7) (其 中 a 设 为 0?) 。 


阵 采 用 D = (OTO)-107 "计算 , WERN Sout =D «Ions o FACES 


E 叶 解 调 方法 时 我 们 使 用 
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1.000 0.709 0.510 —0.298 
1.000 —0.311 0.510 一 0.719 
1.000 —0.311 一 0.510 一 0.719 
1.000 0.709 —0.510 —0.298 
1.000 —0.311 0.510 0.298 
1.000 —0.311 0.510 0.719 
1.000 0.709 -—0.510 0.179 


1.000 0.709 —0.510 0.298 
(3) 初始 方位 角 a 在 为 0"~0.1° 范 围 内 变化 , 根据 步骤 (1) 和 (2) 定量 计算 a 与 解 调 结果 相 
对 误差 的 大 小 关系 。 如 图 3 所 示 。 
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图 4 旋转 波 片 的 初始 方位 角 误差 对 解 调 结果 的 影响 。 己 知 的 输入 偏振 信号 为 Sm = (1 0.02.0.06,0.2)7。 调 
制 方式 详 见 正文 描述 。 初 始 方位 角度 误差 (om) 变化 范围 为 0*~ 0.1*。 我 们 对 观测 信号 分 别 采 用 传 里 叶 分 析 
("x "点 ) 和 解 调 矩 阵 (虚线 ) 两 种 方法 进行 解 调 。 从 上 到 下 分 别 展示 了 两 种 方法 的 解 调 出 的 Q，U，Y 信号 与 

输入 值 之 间 的 相对 误差 。 


Fig. 4 Relative errors of the Stokes parameters Q,U,V caused by the home-position error of the 


rotating waveplate. The input Stokes vectors are set to be Sin = (1,0.02,0.06,0.2)7 . The modulation 


process can be seen in the text. The results derived by the Fourier demodulation are shown by 
the symbols of "x", while the results derived by using the demodulation matrix are plotted by 


the dashed line. 


图 4 展示 了 采用 两 种 解 调 方法 得 到 的 解 调 结果 相对 误差 与 初始 方位 角 的 关系 。 虚 线 代 表 
解 调 矩 阵 方法 的 结果 , 'x' 号 代表 傅 里 叶 分 析 的 结果 。 首先, 无 论 用 哪 种 方法 , 随 着 a 的 增加 ， 
Q@、U ,TY 的 相对 误差 均 明显 增加 。 对 于 线 侦 振 信 号 来 说 ， 两 种 方法 给 出 的 结果 基本 一 致 ( 残 
差 非常 小 ， 在 10- 量 级 以 下 ) 。 但 是 对 圆 侦 振 信号 来 说 ， 傅 里 时 分 析 结 果 呈 现 出 较 大 的 相对 
误差 大约 是 用 解 调和 矩阵 得 到 的 2 倍 左 右 。 和 针对 这 个 情况 我 们 做 了 进一步 分 析 ， 发 现 可 以 遂 
过 增加 一 个 调制 周期 内 的 探测 器 的 采集 频率 (例如 由 8 次 采集 增加 为 16 次 )， 但 是 这 种 方法 在 
实际 工作 中 并 不 太 适 用 ， 我 们 会 在 结果 部 分 进行 讨论 。 其 次 ， 图 中 Q,，U 的 相对 误差 有 一 定 
的 差别 , 说 明 相 对 误差 的 大 小 与 信号 本 身 的 强 弱 有 关 。 当 我 们 进一步 将 输入 的 0，U,，V 信 号 
设 定 为 一 样 大 小 的 时 候 , 他 们 的 相对 误差 则 基本 相同 , 说 明 初 始 方 位 角 误差 对 于 三 者 的 影响 
ve BU. 因此， 对 于 比较 弱 的 10? 量 级 的 偏振 信号 来 说 ,例如 此 处 的 @ 信 号 (Q = 0.02) ， 两 
种 解 调 方 法 都 给 出 一 致 的 要 求 ， 即 ， 当 相对 误差 控制 在 103 量 级 以 内 时 ， 初 始 方 位 角 的 误差 
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要 求 小 于 0.004°( 约 14”") 。 
2.2 波 片 定位 误差 

波 片 旋转 定位 误差 主要 是 由 带动 其 旋转 的 电机 的 转速 不 稳定 造成 , 转速 不 稳定 会 导致 在 
探测 器 每 次 曝光 过 程 中 ， 波 片 快 轴 的 旋转 角度 不 再 是 一 个 稳定 值 ( 即 每 次 的 A91 不 同 ) 。 经 过 
调研 得 知 ， 目 前 常用 的 高 精度 电机 在 达到 匀速 转动 的 状态 之 后 ,重复 定 位 精度 可 以 达到 10” 
左右 ， 因 此 我 们 针对 这 个 量 级 的 定位 误差 进行 了 蒙特 卡 洛 模 拟 。 我 们 的 分 析 思 路 如 下 : 

(1) 输入 已 知 偏振 信号 5;,, = (1, 0.02,0.06,0.2)7。 模 拟 8 步 连续 式 调 制 ， 探 测 器 单 次 曝光 
过 程 中 波 片 转动 角度 rx/8 土 ag，a 为 0 到 10” 内 的 随机 数 。 波 片 延 迟 量 127°%。 根 据 公 式 (6) 生 成 
所 观测 到 的 强度 轮廓 Jops。 


(D 采用 理想 解 调和 矩阵 和 傅 里 叶 分 析 方 法 分 别 对 观测 到 的 强度 轮廓 进行 解 调 。 计 算 解 调 
结果 与 输入 值 的 相对 误差 。 
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图 5 波 片 旋转 定位 误差 对 解 调 结果 的 影响 。 输 入 偏振 信号 与 图 3 相同 。 调 制 方式 详 见 正文 描述 ， 调 制 过 程 
中 波 片 旋转 定位 误差 范围 土 10”"， 随 机 模拟 100 次 独立 调制 过 程 。 分 别 采 用 傅 里 叶 分 析 (" x "点 ) 和 解 调和 所 
PE (虚线 ) 两 种 方法 解 调 并 给 出 解 调 结果 与 输入 值 之 间 的 相对 误差 。 


Fig. 5 Relative errors of the Stokes parameters Q, U, V caused by the angle accuracy of the rotation 


waveplate. Input Stokes parameters are the same as the one in Fig. 3. 100 times modulation are 


operated with the angle accuracy within a range of +10” .The symbols "x" indicate the relative 
errors of the results inferred by the Fourier demodulation, while the dashed line shows the one 


demodulated by using the demodulation matrix. 


5 展示 了 波 片 定位 误差 对 解 调 结果 相对 误差 的 影响 。 在 模拟 了 100 次 独立 调制 的 结果 
后 , 我 们 发 现 : 两 种 解 调 方法 给 出 的 线 偏振 信号 以 及 圆 偏振 的 相对 误差 非常 接近 。 当 输入 信 
号 是 Sm = (1, 0.02,0.06,0.2)7 HY, Q, U, V 的 相对 误差 的 rms (HAP HIN 3.7X 10°. 1.1X10° 
和 1.8X10”。 我 们 进一步 计算 发 现 当 输入 信号 Q =U =V = 0.02 时 ,三 者 的 相对 误差 的 rms 
值 分 别 是 3.7X10 、3.5X10 ”和 2.9X10”。 这 不 仅 说 明 相对 误差 的 大 小 与 信号 的 强 弱 有 
关 ， 也 说 明定 位 误差 对 圆 偏 振 信 号 造成 的 影响 略为 小 于 对 线 偏振 信号 的 影响 。 最 后 ， 对 10” 
量 级 Q 信号 来 说 ，10” 的 波 片 定位 精度 所 产生 的 相对 误差 也 在 10 的 量 级 左右 。 


2. 3 时 间 差 占 比 引起 的 测量 误差 
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图 6 波 片 调制 周期 与 探测 器 采集 频率 的 时 间 差异 。 假 设 波 片 完成 一 个 调制 周期 需要 n xT, E TZAR 
测 器 实际 曝光 时 间 为 Toxj， 两 者 差异 为 A。 


Fig. 6 Difference between the time required for the waveplate to complete one rotation slot during 


one modulation cycle and the time taken by the detector to collect the photons. 

假设 波 片 无 定位 误差 ， 并 在 7 时 间 间 隔 内 匀速 旋转 过 一 个 特定 角度 (180"/m) ， 然 而 ， 期 

间 探 测 器 的 实际 曝光 时 长 只 有 Tsw,， 两 者 存在 一 个 时 间 差 A (Ti — Tw 站， 探测 器 并 没有 采集 

A 时 间 内 的 光 强 , 故而 产生 了 误差 。 我 们 将 两 者 的 时 间 差 A 与 时 间 T 的 比值 称 为 时 间 差 占 比 ( 即 

A/T) 。 此 时 探测 器 单 次 嘿 光 时 间 内 波 片 实际 转动 过 的 角度 为 (1 一人) * 180/n。 在 图 7 中 ,我 

们 模拟 了 A/T 从 0% 增 加 到 2% 时 ， 解 调 结果 与 输入 值 之 间 的 相对 误差 情况 。 在 本 节 中 ， 输 入 
的 偏振 信号 为 Sm = (1 0.02.0.02.0.02)7， 即 线 偏振 与 圆 偏振 信号 强 弱 相等 。 
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Al 7 时 间 差 占 比 对 解 调 结果 的 影响 。 输 入 信和 号 大 小 与 图 3、4 中 的 相同 。 时 间 差 占 比 由 0 增加 至 2%, fie 
则 结果 分 别 采 用 傅 里 叶 分 析 (" x "号 ) 和 解 调 矩阵 (" + "号 ) 获得 。 上 、 中 子 图 中 的 虚线 代表 10°, PFA 
的 虚线 代表 10°. 


Fig. 7 Effect of time difference between the waveplate rotation slot time and the detector frame 


Ss 


time on the modulation. The x-axis means the time difference as a percent of the waveplate rotation 
slot time. "x" represents the results retrieved from the Fourier demodulation, while the "+" 
indicates the one using the demodulation matrix. Dotted lines in the upper and middle panels indicate 


the errors of 10°, while the one in the bottom indicate the error of 10° 


由 图 7 我 们 发 现 : 两 种 解 调 方法 得 到 的 结果 基本 一 致 ， 当 A/T 逐 渐 增 大 时 ， 解 调 结果 的 
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相对 误差 都 会 随 之 增 大 , 两 种 结果 的 明显 区 别 仍然 是 针对 圆 偏振 信号 的 解 调 ， 使 用 解 调和 矩阵 
得 到 的 相对 误差 会 比较 小 。 同 时 我 们 可 以 看 到 ， 相 比 于 圆 偏振 信号 ， 时 间 差 占 比 的 变化 对 线 
偏振 信号 造成 的 影响 更 为 明显 。 从 U 信 号 分 析 发 现 , FA RE EPS HITE 103 的 量 级 时 , A/T 
要 小 于 1%. 我 们 同时 也 分 析 了 Si = (1 0.02,0.02,0.2)7 的 情况 ， 即 了 信号 增强 了 10 倍 ， 依 然 
发 现时 间 差 占 比 的 变化 对 圆 偏 振 信 号 的 影响 也 非常 弱 , 与 图 7 的 结果 非常 类 似 , 因此 没有 展 
示 。 这 说 明 时 间 差 主要 引起 的 是 线 偏振 信号 之 间 的 串扰 。 最 后 这 里 需要 指出 的 一 点 是 ， 本 工 
作 的 模拟 只 是 针对 1 个 调制 周期 ,如 果 进 行 多 个 周期 的 调制 ， 那 么 微小 的 A 会 积 少 成 多 ， 误 
差 会 逐步 放大 。 


上 


3 讨论 与 结论 


本 文 给 出 在 连续 式 旋转 波 片 调制 模式 下 的 传 里 叶 分 析 的 推导 过 程 和 解 调 公式 (公式 7)， 
并 利用 理论 偶 振 轮廓 对 公式 的 准确 性 进行 了 验证 。 这 里 我 们 需要 着 重 指出 公式 (7) 中 对 V 信 
号 的 解 调 ,利用 虚 部 比 实 部 能 得 到 更 加 准确 的 结果 。 接 着 ,我 们 分 别 采用 了 傅 里 叶 分 析 和 经 


= 


E LFS EE RP I Dr E, 利用 理论 模拟 定量 计算 了 波 片 的 初始 方位 角 误差 、 旋 转 定位 误 
差 以 及 探测 器 曝光 过 程 中 时 间 差 占 比 三 个 因素 对 解 调 结果 准确 性 造成 的 误差 , 即 对 偏振 测量 
© 精度 的 影响 。 这 里 我 们 逐一 对 主要 结论 进行 讨论 和 总 结 ; 

O (1) 波 片 初始 方位 角 误差 的 影响 : 对 于 线 偏振 信号 来 说 , 采用 傅 里 叶 分 析 和 解 调 矩 阵 可 
S 以 得 到 比较 一 致 的 结果 。 当 测量 的 偏振 信号 在 102 量 级 时 (这 在 太阳 物理 实测 时 是 常见 的 ) ， 
© 若 要 求 结果 的 相对 误差 在 103 量 级 以 内 ， 那 么 初始 方位 角 的 误差 要 求 小 于 0.004e ( 即 14” 左 
c Ai) 。 但 是 ,我 们 发 现 对 于 圆 偏振 信号 来 说 ， 传 里 叶 分 析 的 方法 得 到 的 相对 误差 比较 大 ， 从 
© 这 一 点 来 看 ,采用 解 调 箱 阵 方法 可 以 克服 初始 方位 角 误差 带 来 的 影响 。 B BRAT SEIN AI 
0 方位 角 误 差 对 圆 偏振 和 线 偏振 产生 的 影响 是 一 致 , 是 相对 误差 的 大 小 与 偏振 信号 本 身 的 强 弱 
= 月 关 。 这 里 我 们 对 比 了 Liang 等 人 (2019) 的 结果 中。 他 们 的 误差 分 析 是 基于 步 进 式 施 转 波 片 
N 调制 及 简单 的 全 里 叶 分 析 方 法 。 他 们 的 结果 (我 们 进行 了 同样 的 模拟 并 肯定 了 他 们 的 结果 


(€ 
w 
a 


V u 


EREA AFA RREN TBA de r AA RZE SLP i j 影 响 。 他 们 基于 线 俩 振 信 号 的 分 析 
结果 则 与 本 文 一 致 。 因 此 我 们 认为 ， 相对 于 步 进 式 调 制 ， 连 续 式 调制 对 初始 方位 角 的 定位 精 
度 要 求 更 高 。 


(2) 波 片 旋转 定位 误差 的 影响 : 无 论 对 于 线 偏振 还 是 圆 偏 振 信和 号， 两 种 解 调 方法 结果 基 
本 一 致 。 两 个 结果 均 显 示 旋 转 定位 误差 对 线 偏振 信号 造成 的 影响 更 为 明显 ， 且 相对 误差 的 大 
小 与 信号 的 强 弱 有 关 。 定 量 分 析 发 现 ， 当 波 片 的 重复 定位 精度 在 10” 左 右 时 (符合 常见 高 精 
度 电机 的 性 能 指标 )， 对 107? 量 级 的 弱 偏 振 信号 来 说 , 测量 误差 也 能 在 103 的 量 级 左右 。Liang 
等 人 (2019)"“ 给 出 在 步 进 式 调制 时 ， 当 要 求 相对 测量 误差 在 10-3 的 量 级 左右 时 , 波 片 的 定位 
精度 可 以 放宽 到 0.1?。 因 此 我 们 认为 在 相同 的 测量 误差 要 求 的 前 提 下 ， 连 续 式 调制 比 步 进 式 
调制 模式 对 波 片 旋转 的 定位 精度 要 求 更 高 。 


(3) 探测 器 采集 过 程 中 时 间 差 占 比 的 影响 : 这 个 影响 仅仅 存在 于 连续 式 调制 模式 下 。 在 
理想 情况 下 ， 我 们 要 求 7 = Terp ,而 实际 测量 情况 下 ， 由 于 图 像 存 储 或 数据 转移 需要 一 定 的 
时 间 A， 这 会 造成 调制 方式 偏离 理想 情况 ， 导 致 解 调 结果 出 现 误差 。A 一 般 在 hs 量 级 ， 相 对 
PAIRS RRI F ms ARE CRT TLR Oe 这 个 值 非常 小 。 而 且 我 们 发 现 ， 时 间 差 造 
成 的 影响 主要 是 线 偏振 之 间 的 相互 串扰 , AT 在 1% 左 右 时 ,对 圆 偏振 信号 的 相对 误差 的 影响 
接近 10" 量 级 。 另 外 ， 在 连续 式 采集 模式 下 ， 除 了 时 间 差 以 外 ， 我 们 还 应 关注 探测 器 采集 帧 
频 与 波 片 旋转 速度 之 间 的 匹配 问题 , 利用 探测 器 的 采集 帧 频 对 波 上 方 的 旋转 速度 下 限 提 出 要 求 。 


(4) 最 后 值得 指出 的 是 , 虽然 没有 定量 计算 , 但 我 们 发 现 ， 增加 一 个 调制 周期 内 探测 器 
的 采集 频率 后 ， 无 论 使 用 哪 种 解 调 方法 ， 侦 振 测量 精度 都 会 有 所 提高 ， 但 是 增加 采集 频率 会 
造成 数据 量 的 增加 以 及 调制 频率 或 者 时 间 分 辩 率 的 降低 " ， 因 此 综合 考虑 各 种 因素 后 ， 我 
们 在 本 文 的 工作 中 仍 采用 常见 的 一 个 周期 采集 8 帧 的 模式 。 
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Fourier Demodulation of the Continuously Rotating 
Waveplate Modulation and the Error Analysis 


Lin Zheyu'”, Xu Zhi 
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Abstract: Continuously rotating waveplates as the polarization modulator is an 


important candidate polarimeter configuration for many solar telescopes to measure the 


solar magnetic field. It is widely adopted in many ground-based or space-based solar 


telescopes for the solar magnetic field measurement. In this paper, we present the 


derivation and formula of the Fourier demodulation in the context of continuously rotating 


waveplates. We further testify the formula correctness by using the synthetic Fel Stokes 


profiles, which are produced by the RH radiative transfer code with a given atmosphere 


model. In addition, we calculate the effect of the home-position errors, the positioning 
errors of the waveplate rotation and the time difference between the waveplate slot time 
and the detector frame time on the demodulated results, by using two kinds of 
demodulation methods, i.e., the demodulation matrix and the Fourier analysis. The main 
findings are: (1) The new Fourier analysis can provide more accurate results than the 
simple one adopted by the step-wise modulation. Similar results are obtained when we do 
the relative error estimation using these two demodulation methods. (2) Considering the 
relative errors due to the home-position angle error, we find the demodulated linear 
polarization are almost the same using these two methods. However, in the case of 
circular polarization, the demodulation based on the demodulation matrix is more reliable. 
Whatever which method is applied, it is shown that the home-position angle error has the 
same effect on both the circular and linear polarization. The relative error of 10-3 requires 
the home-position angle accuracy within ten arc seconds. (3) Calculating the relative error 
caused by the position error of the waveplate rotation, we find the demodulation results 
are almost equal using these two methods. Both of them show that the position error has 
more influence on the linear polarization. Furthermore, the requirements of the position 
accuracy are much higher (ten of arc seconds) than that in the case of step-wise 
modulation (0.1 degree) in order to achieve the relative error of 10-3. (4) It is crucial to 
precisely match the detector frame time and the waveplate rotation slot time. We 
investigate the relative error due to the difference between them and find that the time 
difference (A) can cause cross-talk between the linear polarization signals. It has more 
prominent effect on the linear polarization than on the circular polarization. The time 
difference as a percent of the waveplate slot time (A/T) is required to be smaller than 1% 


if the relative error is smaller than 10°3 for the linear polarization. 
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